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ABSTRACT 

(Unclassified) 

The  generation  of  a  high— intensity  acoustic  field  which  is  a  faithful  replica, 
or  at  least  a  known  function,  of  the  input  signal  is  a  prime  requirement  in  modern 
active  sonar  systems.  Most  sonar  signals  in  use  are  modulated  sine  waves. 
Therefore,  transducers  must  be  capable  of  efficiently  transforming  this  type  of 
signal.  It  is  evident  that  design  criteria  and  test  procedures  based  on  steady-state 
sine-wave  signals  are  not  adequate.  The  potential  gain  in  sonar  performance  to  be 
realized  by  maximizing  signal  power  without  distortion  merits  increased  effort 
addressed  to  this  problem. 

This  report  is  intended  to  point  out  some  of  the  problems  associated  with  the 
generation  of  maximum-energy  acoustic  fields  by  large  multielement  sources  and  to 
suggest  techniques  of  analysis  and  testing.  The  use  of  electrical  analog  analysis  of 
power-limiting  factors  in  transducer  elements  appears  to  be  an  attractive  method  if 
a  means  can  be  developed  for  estimating  the  range  of  acoustic-loading  values 
experienced  by  e  lements  in  an  active  array . 

A  description  of  calibration  procedures  used  and  results  obtained  in  tests  of 
the  ARTEMIS  acoustic  source  are  presented  to  serve  as  an  example  of  proposed  test 
methods . 


PROBLEM  STATUS 

This  is  an  interim  report  on  one  phase  of  this  project.  Work  is  continuing . 
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INTRODUCTION 

Most  acoustic  generators  of  practical  interest  in  sonar  employ  modulated  sine- 
wave  signals  to  convey  Information.  The  simplest  form  of  modulation  Is  widely  used 
in  echo  ranging  and  consists  of  a  sinusoidal  signal  which  is  gated  on  periodically  to 
produce  pulses  of  energy  having  typical  time  durations  corresponding  to  a  large 
number  of  cycles  of  the  modulated  wave.  The  probability  of  detecting  a  signal  in  the 
background  noise  increases  with  increased  signal  energy.  Signal  energy  con  be 
increased  by  increasing  signal  power  or  duration,  or  both.  When  a  target  echo 
contains  sufficient  doppler  shift  to  permit  discrimination  between  signal  and  rever¬ 
beration,  the  simple  gated  sine-wave  signal  is  well  suited  for  detection. 

Increasing  pulse  duration  and  power  enhances  performance.  Accurate  location  of 
the  target  in  range,  however,  sets  a  limit  on  the  pulse  duration  depending  on  the 
resolution  requirements.  Beyond  that  point,  signal  energy  can  be  increased  only  by 
Increasing  power  or  by  Increased  focusing  of  the  energy  in  the  direction  of  the 
target.  The  use  of  nonperiodic  signals,  together  with  coherent  processing,  eliminates 
the  dependence  of  bandwidth  on  pulse  length  and  therefore  permits  the  use  of  long- 
duration  signals  with  acceptable  range  resolution.  Such  signals  are  also  effective 
against  low  doppler  targets.  Signal  energy  cannot  be  Indefinitely  increased  by 
Increasing  signal  duration  due  to  the  temporal -stability  limitation  of  the  propagation 
path.  Therefore,  with  either  gated  sine-wave  or  the  nonperiodic  signals,  it  is 
desirable  to  maximize  signal  power.  In  systems  employing  coherent  processing,  this 
must  be  accomplished  with  a  minimum  of  phase  and  amplitude  distortion  over  the 
frequency  band  containing  the  signal . 

As  the  power  Input  to  a  transducer  Is  Increased,  all  factors  in  the  transducer 
element,  such  as  potentials,  currents,  temperatures,  pressures,  and  mechanical 
stress,  approach  their  limiting  values,  which,  If  exceeded,  will  either  cause  failure 
of  the  element  or  cause  It  to  operate  in  an  unacceptably  nonlinear  fashion.  That 
factor  which  first  reaches  its  limit  defines  the  maximum  allowable  power  input  to 
the  transducer.  If  operating  conditions,  such  as  acoustic  loading,  ambient  pressure 
or  temperature,  or  type  of  signal  used,  are  changed,  there  may  be  a  change  in  the 
power-limiting  factor.  All  the  power-limiting  factors,  with  the  exception  of 
certain  thermal  conditions,  are  sensitive  to  peak  rather  than  average  values.  In 
general,  the  transient  response  of  a  single  element  or  an  element  operating  in  an 
array  of  elements  is  such  that  It  Is  difficult  to  determine  the  peak  values  of  power- 
limiting  factors  analytically  and,  In  most  cases,  it  is  equally  difficult  to  measure 
them  experimentally.  In  practice,  the  failure  of  an  element  Is  frequently  the  first 
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Indication  that  a  limiting  value  has  been  reached.  Information  obtained  In  this 
manner,,  hnw*y»r,  net  necessarily  vollu  if  operaring  conditions  have  changed. 

On  the  other  hand,  if  transducer  elements  are  not  operated  at  power  levels  close  to 
their  limit,  system  performance  is  degraded  by  the  use  of  less-than-maximum  avail¬ 
able  signal  power.  Clearly,  a  need  exists  for 

1 .  Increasing  knowledge  of  the  causes  of  power  limitations  in  acoustic  sources, 

2.  Developing  nondestructive  techniques  of  determining  power  limitations  in 
existing  transducers, 

3.  Establishing  design  criteria  for  optimizing  power  capability  of  sources  to  be 
used  with  particular  signal  types,  and 

4.  Specifying  techniques  for  testing  sources  for  phase  and  amplitude  distortion. 

This  report  is  intended  to  call  attention  to  the  importance  of  these  four  areas 
to  optimized  sonar  performance  and  to  report  on  limited  investigations  which  lead 
to  conclusions  and  recommendations  on  certain  aspects  of  the  problem  and  to  point 
out  examples  of  problem  areas  In  other  facets. 

EFFECT  OF  ACOUSTIC  LOADING  ON  POWER  LIMITATION 

General 

An  acoustic  projector  having  dimensions  which  are  appreciable  compared  to 
the  sound  wavelength  In  water  cannot,  ip  general,  experience  a  uniform  acoustic 
load  over  its  radiating  surface  due  to  acoustic  interactions.  Since  all  large  projec¬ 
tors  are  composed  of  a  large  number  of  individual  transducer  elements,  there  Is  a 
variation  In  loading  from  element  to  element.  In  the  steady-state  condition,  the 
loading  pattern  is  fixed  and  can  be  calculated,  although  it  might  be  necessary  to 
make  various  simplifying  assumptions  In  order  to  restrict  the  complexity  of  the 
calculation  to  a  form  compatible  with  available  computers.  Figure  lz  illustrates 
the  steady-state  loading  pattern  calculated  for  one  quadrant  of  a  2-1/2^  by  4\ 
rectangular  plane  array  In  which  the  simplifying  condition  of  a  uniform  velocity 


1 .  R.  V.  Baler,  "Theoretical  Interaction  Computations  for  Transducer  Arrays, 
Including  the  Effects  of  Several  Different  Types  of  Electrical  Terminal  Connec¬ 
tions",  NRL  Conf  Report  6314,  Oct  1965 

2.  R.  H.  Ferris,  R,  V.  Baler,  "Acoustic  Interaction  In  Multielement  Sonar  Arrays" , 
Proceedings  of  the  Seventh  Navy  Science  Symposium,  ONR-16,  vol,  2,  May 
1963  (Conf) 
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over  the  face  of  the  array  was  assumed.  It  can  be  seen  that  very  small  elements  In 
such  an  array  would  experience  wide  variations  in  loading,  depending  on  their 
position  in  the  array.  Larger  elements  would  tend  to  average  the  loading  on  their 
individual  radiating  faces  and  thus  would  be  subjected  to  a  lessor  variation.  It  Is 
possible,  however,  that  the  larger  elements  could  be  adversely  affected  by  unsym- 
metric  loading.  In  contrast  to  the  steady-state  problem,  a  complete  analytical 
solution  for  the  conditions  encountered  with  modulated  signals  is,  in  the  presont 
state  of  the  art,  unthinkable.  Fortunately,  a  complete  solution  is  not  required  in 
order  to  determine  the  power-limiting  effects  of  the  loading  variations.  It  is  suf¬ 
ficient  to  know  the  range  of  possible  complex  loadings.  When  the  behavior  of  an 
element  is  examined  over  this  range,  the  worst  case  can  be  assumed  to  be  /he 
limiting  condition,  since  the  allowable  power  Input  to  the  array  1$  limited  by  the 
worst-case  element.  Although  a  technique  of  estimating  the  range  of  loading 
encountered  with  modulated  signals  is  not  in  hand,  it  appears  to  be  a  more  tractable 
problem  than  a  complete  analytical  solution. 

A  possible  approach  to  investigating  the  effects  of  loading  on  power-handling 
capability  is  to  Instrument  a  sample  of  elements  to  record  the  instantaneous  response 
of  each  of  the  factors  which  might  limit  power  while  the  full  array  Is  being  driven 
with  the  particular  modulated  signals  to  be  used  In  the  system  of  which  the  projector 
is  a  part.  This  approach  has  some  obvious  disadvantages  In  that  the  required  Instru¬ 
mentation  would  generally  require  a  major  modification  of  the  elements  and  tests 
must  be  made  on  the  completed  source.  A  second  approach  would  be  to  Instrument 
a  single  element  on  which  the  dynamic  load  can  be  varied  over  the  required  range 
by  some  technique  such  as  the  standing -wave  tube  or  dynamically  driven  load 
elements.  A  third  approach  which  is  applicable  to  some  transducers  involves  the 
use  of  an  electrical  analog  circuit.  If  such  a  circuit  can  be  constructed  to  accurately 
represent  an  element,  measurement  of  the  power-limiting  factors  is  greatly  simplified 
and  loading  Is  completely  and  readily  controllable.  An  example  of  this  technique 
will  now  be  shown. 

Analog  Technique 

Figure  2  Illustrates  the  basic  configuration  of  an  Illustrative  ceramic  element. 

It  consists  of  a  piezoelectric  ceramic  stack  having  a  reaction  mass  bonded  to  one 
end  and  a  radiating  mass  bonded  to  the  other  end.  A  third  mass,  the  watertight 
housing,  is  joined  to  the  radiating  mass  with  a  pressure-release  gasket.  The  ceramic 
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stack  acts  as  a  spring  between  the  reacting  and  radiating  masses,  and  the  pressure- 
release  gasket  acts  as  a  spring  between  the  radiating  mass  and  the  housing.  Figure  3 
illustrates  the  equivalent  circuit  of  this  element  being  driven  from  an  amplifier  hav¬ 
ing  an  internal  impedance  R.  The  analog  circuit  was  assembled  and  measurements 
made  of  various  possible  power-limiting  factors.  In  this  circuit  the  analog  voltage 
is  equivalent  to  force  and  current  to  velocity.  To  exemplify  the  analysis  technique, 
several  possible  power-limiting  factors  were  selected  and  their  amplitudes  recorded 
as  the  analog  circuit  was  driven  by  continuous-wave  and  modulated  waveforms. 

The  specific  factors  selected  were 

1 .  The  mechanical  stress  in  the  ceramic,  given  by 


where 

V  Is  the  voltage  across  capacitor  c. , 

C1  1 

A  is  the  area  of  the  ceramic  material,  and 
c 

9  Is  the  electromechanical  turns  ratio. 

2.  The  mechanical  stress  in  the  pressure-release  material  between  the  radiating 
and  housing  masses,  given  by 

where 

V  Is  the  voltage  across  capacitor  c,  , 

c3  J 

A^  is  the  area  of  the  pressure-release  material,  and 

9  is  the  electromechanical  turns  ratio. 
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3.  The  potential  gradient  In  the  ceramic,  given  by 
V 

c 

F  =  _ “ 

8  d  ' 

where 

V  is  the  voltage  across  capacitor  c  ,  and 
o 

d  is  the  thickness  of  one  ceramic  disk. 


For  illustrative  purposes,  maximum  allowable  values  for  each  of  these  three 
factors  were  arbitrarily  selected  and  corresponding  peak  values  of  electrical  poten¬ 
tials  across  the  appropriate  components  in  the  analog  circuit  were  computed  from  the 
given  equations.  Each  of  the  three  potentials  was  held  constant  in  turn  at  the  com¬ 
puted  values  as  the  circuit  was  driven  by  steady-state  continuous-wave  signals  which 
were  stepped  in  frequency  across  the  operating  band.  During  this  procedure  the 
voltage  input  was  recorded.  The  assumed  limiting  values  were  1 25-pounds-per- 
square-tnch  stress  in  the  ceramic,  25-pounds-per-square-lnch  stress  In  the  pressure- 
release  material,  and  a  maximum  potential  E  across  the  ceramic  stack  of  1000  volts. 


This  process  was  repeated  for  various  equivalent  values  of  acoustic  load,  as 
given  in  Table  1 . 


TABLE  1 


Real  Part/ocA 

Reactive  Part/pcA 

0.0 

-0.5 

0.0 

0.0 

0.0 

+0.5 

0.1 

-0.5 

0.1 

0.0 

0.1 

+0.5 

0.5 

-0.5 

0.5 

0.0 

0.5 

+0.5 

1.0 

-0.5 

1.0 

0.0 

1.0 

+0.5 
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A  positive  reactive  load  represents  mass  loading  and  a  negative  sign  represents 
stiffness  loading  ■ 

The  results  are  plotted  in  Fig.  4,  where  voltage  input  is  plotted  as  a  function 
of  frequency.  These  curves  illustrate  the  maximum  allowable  steady-state  input 
voltage  as  determined  by  the  worst-case  loading  condition  of  the  three  factors 
investigated.  The  curves  in  Fig.  4  represent  the  power  limitation  imposed  by  the 
factors'  stress  in  the  ceramic,  stress  in  the  band,  and  potential  gradient  in  the 
ceramic  .  Although  the  power  output  of  an  element  depends  upon  the  acoustic  load 
seen  by  that  particular  element  as  well  as  the  driving  voltage,  the  total  power 
output  of  an  array  of  parallel-connected  elements  having  a  fixed  loading  distribu¬ 
tion  is  limited  by  the  maximum  allowable  voltage  drive  to  the  worst-case  element. 

It  can  be  observed  that,  for  the  assumed  limiting  values  of  the  three  factors,  the 
power  output  is  not  always  limited  by  the  same  factor.  If,  as  in  this  example,  it 
was  determined  that  the  steady-state  loading  at  various  points  in  an  array  of  these 
elements  varied  over  the  limits  of  0.0  to  1 .0  for  the  real  part  and  -0.5  to  0.5  for 
the  reactive  part,  the  power  Input  would  be  limited  by  potential  gradient  between 
4.0  and  4,43  kilohertz,  by  stress  in  the  ceramic  between  4.43  and  4.72  kilohertz, 
and  by  stress  In  the  pressure  release  for  higher  frequencies. 

A  comparable  analysis  for  modulated  waveforms  is  complicated  by  the  fact  that 
with  modulated  signals  the  interaction  field  seen  by  an  element  Is  nonperiodic.  A 
possible  approach  is  to  estimate  the  range  over  which  the  load  can  vary  both  in 
time  and  position  and  provide  a  dynamic  load  for  the  analog  circuit  which  varies 
in  a  random  fashion  within  the  limits  of  this  range.  The  modulated  signal  Input  to 
the  analog  circuit  could  then  be  increased  in  value  until  the  peak  value  of  the 
factor  being  investigated  reached  its  limiting  value  at  least  once.  It  is  possible 
that  the  value  of  the  maximum  allowable  input  voltage  obtained  in  this  statistical 
manner  would  be  somewhat  lower  than  that  required  in  a  real  array.  On  the  other 
hand,  an  approximate  analysis  can  be  made  utilizing  a  selection  of  fixed  loads,  as 
was  done  in  the  steady-state  case.  The  approximation  could  be  expected  to  degrade 
as  the  signal  bandwidth  is  increased.  The  method  employing  dynamically  varying 
loads  has  not  been  investigated  experimentally.  However,  measurements  of  maxi¬ 
mum  input  voltage  were  obtained  for  modulated  signals  employing  a  selection  of 
fixed  loads  to  compare  with  the  steady-stale  measurements  performed  with  the  analog 
circuit  shown  in  Fig.  3.  Pulsed  continuous-wave  and  phase-reversal  modulated 
pseudorandom  sequences  in  which  the  modulation  code  was  generated  by  a  linear 
shift  register  were  used.  Ceramic  strain  was  assumed  to  be  the  limiting  factor. 
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The  input  level  was  increased  until  the  peak  value  of  the  ceramic  strain  reached 
the  limiting  value,  at  which  time  the  root-mean-square  voltage  input  was  recorded. 
The  results  are  listed  in  Table  2.  The  carrier  frequency  in  each  case  was  five  kilo¬ 
hertz.  The  first  column  lists  tho  type  of  modulation  used,  the  second  column  the 
resistive  acoustic  load,  and  the  third  column  the  reduction  in  allowable  input 
voltage  when  compared  with  continuous-wave  operation  under  the  same  conditions. 

TABLE  2 


Type  of  Modulation 

Acoustic  Load 

Reduction  in 

Allowable  Input  Level 
Ref.  to  Steady  State 

Pulsed  continuous  wave 

1 .0  pcA 

3.0  db 

Pulsed  continuous  wave 

0. 1  pcA 

5.0  db 

Phase-modulated  PRS 

1 .0  pcA 

6.0  db 

Phase-modulated  PRS 

0.1  pcA 

1 1 .8  db 

Results  of  Measurements  on  the  ARTEMIS  Source 
3 

During  initial  tests  of  the  partially  completed  ARTEMIS  acoustic  source  in 
May  1961,  more  than  40  percent  of  the  transducer  elements  failed  when  the  array 
was  driven  at  a  power  level  below  the  design  value.  At  the  time  of  this  test,  only 
one-tenth  of  the  elements  had  been  installed,  resulting  in  an  array  approximately 
one  wavelength  square.  In  the  Investigation  of  the  cause  of  the  element  failures, 
accelerometers  were  employed  to  obtain  information  on  the  displacement  amplitudes 
of  the  radiating  faces  of  the  individual  elements.  A  large  variation  in  displacement 
amplitudes  was  observed  with  many  of  the  elements  having  displacements  far  in 
excess  of  predicted  values  based  on  radiated  power  of  the  array.  The  ARTEMIS 
elements  are  of  the  variable-reluctance  type  having  two  masses  and  connecting 
springs.  The  element  radiating  face  is  a  square  plane  with  the  length  of  the  sides 
approximately  1/12  wavelength  In  water.  Examination  of  the  damaged  elements 
revealed  that  the  transducer  springs  had  broken  due  to  excessive  deflection. 
Computer  analysis,  which  was  performed  at  a  later  date,  disclosed  that,  with 
elements  of  this  small  size,  a  large  variation  in  acoustic  loading  among  the  elements 

3.  A.  T.  McClinton,  R.  H.  Ferris,  W.  A.  Herrington;  " Project  ARTEMIS  High 
Power  Acoustic  Source,  Interim  Report  on  Acoustic  Performance";  NRL  Memo¬ 
randum  Report  1205,  Aug  1961  (Conf) 
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could  be  expected  due  to  acoustic  interactions.  The  deleterious  effect  of  inter- 
or  t'or>«^  was  heightened  fey  the  electrical  connection  of  ulemenu  in  which  the 
elements  were  divided  into  groups  of  six  with  a  series  connection  of  elements  within 
each  group  and  a  parallel  connection  of  groups.  Since  the  input  impedance  of  an 
element  increases  when  the  element  becomes  unloaded,  the  series  connection  resulted 
in  increased  power  being  applied  to  unloaded  elements  in  a  group,  thus  further 
increasing  their  deflections.  A  subsequent  modification'1  placed  all  elements  In  a 
parallel  connection,  and  this  greatly  improved  the  velocity  distribution  but  not  to 
the  extent  where  full  design  power  could  be  applied  to  the  array.  A  computer 
analysis  indicated  that  the  real  part  of  the  radiation  loading  in  a  one-wavelength- 
square,  parallel-connected  array  would  vary  in  the  range  from  -0.1  to  +0.9pcA. 

The  computation  was  performed  also  for  the  full-scale  2—  1/2X — by— array .  However, 
in  this  case,  the  computer  and  computation  techniques  utilized  necessitated  assum¬ 
ing  an  equal  velocity  for  all  elements.  The  latter  computation  resulted  in  a  range 
of  values  for  the  real  part  of  the  radiation  loading  of  0.2  to  1 .01  pcA.  All  compu¬ 
tations  were  based  on  steady-state  sinusoidal  drive. 

It  was  evident  from  the  observed  failures  that  the  critical  factor  which  Imposed 
a  limitation  on  operating  power  level  was  the  transducer-element  springs.  No  in¬ 
formation  was  available  as  to  the  actual  deflections  experienced  by  the  springs 
under  various  loading  conditions  and  types  of  input  signals  other  than  that  the  per¬ 
missible  deflections  were  being  exceeded  at  an  unpredlcted  low  level  of  drive  as 
evidenced  by  the  failures.  A  program  was  then  initiated  to  investigate  the  actual 
spring  deflections  by  modifying  a  group  of  elements  by  the  installation  of  accelero¬ 
meters  on  the  two-element  masses  in  appropriate  positions,  such  that  the  spring 
deflections  could  be  computed  from  the  observed  relative  motion  of  the  masses. 
Installation  of  the  accelerometers  required  disassembly  and  reassembly  of  the  elements, 
since  the  springs  and  the  reaction  mass  are  completely  enclosed  within  the  radiating 
mass,  which  serves  as  an  outer  housing .  Tests^  were  made  of  the  completed 
ARTEMIS  source  in  which  these  specially  Instrumented  elements  were  Inserted 
sequentially  in  approximately  150  locations  scattered  over  the  extent  of  the  array. 
Steady-state  measurements  were  performed  at  small  increments  of  frequency  In  the 
range  from  350  to  500  hertz.  The  operating  band  for  the  ARTEMI5  source  is  350  to 
450  hertz.  The  mean  and  maximum  observed  spring  deflections  are  plotted  in  Fig.  5. 

It  can  be  seen  that  the  ratio  of  the  maximum  to  the  mean  is  particularly  high  at  the 
upper  end  of  the  frequency  band  and  at  frequencies  in  the  neighborhood  of  400  hertz. 


4.  R.  H.  Ferris;  "Project  ARTEMIS  High  Power  Acoustic  Source,  Fourth  Interim 
Report  on  Acoustic  Performance",  NRL  Memorandum  Report  1400,  Mar  1963  (Conf) 

5.  R.  H.  Ferris;  "Test  of  Project  ARTEMIS  Acoustic  Source",  NRL  Memorandum 
Report  1648,  Sep  1965  (Conf) 
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The  ratio  at  450  hertz  is  approximately  2-1/2  to  1,  which  results  in  a  reduction  in 
allowable  Input  power  of  about  eight  decibels  relative  to  the  power  which  would 
be  acceptable  if  all  spring  deflections  were  uniform  and  equal  to  the  mean.  Other 
studies  have  shown  that  the  deflection  variation  is  principally  due  to  the  nonuni¬ 
form  loading,  although  a  rotary  motion 3  of  the  inner  mass  accounts  for  some  of  the 
spread.  The  rotary  motion  is  particularly  pronounced  at  frequencies  in  the  neigh¬ 
borhood  of  400  hertz. 

The  maximum  permissible  spring  deflection  is  approximately  0.010  inch  peak 
to  peak.  At  higher  deflections  they  are  subject  to  fatigue  failures.  Applying  the 
spring  deflection  limitation  to  the  data  of  Fig  .  5,  the  data  of  Fig  .  6  were  derived. 
This  figure  illustrates,  as  a  function  of  frequency,  the  relative  permissible  voltage 
input  to  the  amplifiers.  Absolute  values  of  voltage  input  are  not  given,  since  the 
preamplifier  gain  can  be  varied  over  wide  limits.  Finally,  the  maximum  safe  source 
level  can  be  computed  for  steady-state  operation  by  combining  acoustic  calibration 
data  with  the  spring  deflection  data.  This  information  is  plotted  in  Fig.  7.  A 
similar  plot  can  be  made  for  the  maximum  permissible  ppwer  Input  to  the  transducer. 
The  power  Input  data  have  been  supplied  to  the  users  of  the  source  as  a  guide  to 
prevent  damage  to  the  elements.  A  suitable  safety  factor  has  been  included  to 
account  for  error  due  to  the  size  of  the  statistical  sample  from  which  these  data 
were  derived  and  to  allow  for  overshoots,  which  it  will  be  shown  can  occur  when 
the  signal  is  gated  on  or  off. 

The  analog  technique  of  analysis  described  In  an  earlier  section  also  was  applied 
to  the  ARTEMIS  source  element.  The  equivalent  analog  circuit  shown  In  Fig.  8 
was  constructed  and  tested  in  the  laboratory  under  steady-state  conditions.  The 
range  of  values  for  acoustic  loading  was  estimated  to  be  from  0  to  1 .0  pcA  for  the 
real  part  and  from  -0.3  to  +0.3  pcA  for  the  reactive  part.  The  simulated  load  on 
the  analog  circuit  was  ranged  over  these  values  as  the  frequency  was  stepped  over 
the  operating  band.  The  input  voltage  was  maintained  at  values  which  produced  a 
constant  simulated  spring  deflection  of  7.07  mils  peak  to  peak.  This  value  of  spring 
deflection  results  from  applying  o  three-decibel  safety  factor  to  the  limiting  value 
of  ten  mils.  The  value  of  the  Input  voltage  required  to  maintain  the  constant  spring 
deflection  was  recorded.  At  each  frequency,  the  minimum  value  of  input  voltage, 
depending  on  load,  was  selected  and  input  powers  to  the  array  calculated  corres¬ 
ponding  to  the  selected  element  voltages.  This  power  represents  the  predicted 


i,  "See-  footnote  5  on  p.  8." 
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maximum  allowable  input  power  to  the  array.  The  results  are  plotted  in  Fig.  9, 
along  with  the  curve  of  maximum  allowable  Input  power  determined  from  experi¬ 
mental  measurements  with  specially  instrumented  elements  in  the  field.  The  third 
curve  illustrates  the  predicted  maximum  allowable  input  power  for  the  condition 
of  uniform  loading  having  a  value  of  unity  pcA  with  no  reactive  loading.  The  dip 
in  the  experimental  curve  at  400  cycles  per  second  results  from  the  excitation  of  a 
rotary  mode  of  vibration  in  the  element.  This  exemplifies  a  shortcoming  in  the 
analog  circuit  analysis  in  that  the  equivalent  circuit  did  not  include  this  degree  of 
freedom.  The  two  curves  are  In  reasonable  agreement  otherwise  when  it  is  con¬ 
sidered  that  they  ore  both  subject  to  inherent  errors.  The  principal  uncertainty  in 
the  experimental  curve  results  from  the  small  size  of  the  statistical  sample  of 
elements;  i.e.,  approximately  ten  percent  of  the  elements  were  instrumented.  On 
the  other  hand,  the  analog  data  suffer  from  the  lack  of  certainty  concerning  the  range 
and  combination  of  real  and  Imaginary  components  of  acoustic  loading  values  and 
ihe  approximation  of  the  equivalent  circuit  to  the  actual  element.  For  the  analog 
technique  to  have  real  value  in  predicting  available  power  levels,  a  method  must 
be  developed  for  accurately  estimating  the  limiting  bounds  of  acoustic  loading. 

The  maximum  allowable  Input  power  to  the  ARTEMIS  source  as  defined  in  Fig.  9 
refers  to  steady-state  sinusoidal  drive.  This  type  of  operation  has  no  practical 
Interest  in  sonar,  nor  can  it  be  assumed  that  the  same  power  limitations  apply  to 
modulated  signals.  The  spring  deflections  of  20  elements  of  the  ARTEMIS  source 
array  were  monitored  while  the  array  was  driven  with  a  shift-register-generated 
pseudorandom  sequence.  It  was  observed  that  occasional  large  peak  spring  deflec¬ 
tions  occurred  In  each  element.  Measurement  of  these  peak  values  disclosed  that 
the  maximum  safe  Input  power  for  this  type  of  signal  was  less  than  the  maximum 
allowable  steady-state  Input  with  a  sinusoidal  signal .  The  pseudorandom  signal 
(I.e.,  PRS)  employed  was  of  the  phase-reversal  modulated  type  which  has  a  peak- 
to-rms  ratio  only  slightly  greater  than  a  sine  wave.  The  maximum  allowable  PRS 
powers,  as  determined  from  experimental  measurements  on  the  array  and  from  the 
analog  circuit  tests  with  the  assumed  loading  are  plotted  as  points  in  Fig  .9.  It  can 
be  seen  that  the  two  methods  agree  within  one  decibel .  This  is  surprisingly  good 
agreement,  considering  that  only  fixed  loads,  whose  range  of  values  was  derived 
from  steady-state  computations,  were  used.  Further  studies  concerning  the  nature 
of  transient  loading  effects  are  needed  before  conclusions  can  be  drawn  regarding 
the  adequacy  of  this  technique. 
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PHASE  AND  AMPLITUDE  DISTORTION 
General 

The  effectiveness  of  coherent  processing  techniques  depends  upon  having 
complete  knowledge  of  the  transmitted  signal .  In  addition  to  the  normal  calibra¬ 
tion  procedures,  the  complex  transfer  function  should  be  measured  for  any  acoustic 
source  which  Is  to  be  used  In  a  coherent  system.  The  measurements  should  be  made 
on  the  entire  source  system,  including  filters,  preamplifier,  power  amplifier,  load¬ 
matching  components,  and  transducer.  Ideally,  the  entire  system  should  be  tested 
together  under  operating  conditions.  By  repeating  the  measurements  at  several 
power  levels  from  less  than  one  percent  of  full  power  to  full  power,  the  linearity 
can  be  determined.  If  the  system  is  linear,  sine-wave  measurements  at  small 
increments  of  frequency  throughout  the  operating  band  will  suffice  to  define  the 
system  characteristics  for  all  signal  types.  On  the  other  hand,  If  significant  non¬ 
linearity  Is  observed,  the  transfer  function  must  be  obtained  by  spectral  analysis 
with  the  types  of  signals  of  Intended  use. 

The  measurement  of  source  transfer  function  requires  that  a  monitoring  hydro¬ 
phone  be  held  In  a  fixed  position  relative  to  the  acoustic  source  and  the  position 
be  accurately  known.  At  high  frequencies,  the  required  degree  of  stability  Is 
difficult  to  achieve.  For  example,  a  system  operating  at  a  frequency  of  five  kilo¬ 
hertz  would  require  that  the  hydrophone  maintain  Its  position  constant  relative  to 
the  source  to  within  one-sixth  of  an  inch  In  range  in  order  to  maintain  a  phase 
stability  of  five  degrees.  The  hydrophone  distance  from  the  source  must  be  suffi¬ 
cient  to  obtain  an  approximation  of  far-fteld  characteristics  to  the  desired  degree 
of  accuracy . 

The  transfer  function  of  the  ARTEMIS  acoustic  source  was  measured  at  sea  during 
November  1965.  During  these  tests  the  source  ship,  the  USNS  MISSION  CAPIS¬ 
TRANO  (T-AG  162),  was  lying  to  In  the  sheltered  waters  of  Northwest  Providence 
Channel  with  the  transducer  array  submerged  to  a  depth  of  600  feet.  Three  hydro¬ 
phones  were  mounted  on  a  spar  which  was  affixed  to  the  end  of  a  190-foot  boom. 

The  boom  was  pivoted  at  the  base  of  the  transducer  array  in  such  a  manner  that  the 
hydrophones  could  be  held  rigidly  In  front  of  the  array  In  several  fixed  positions  In 
a  vertical  plane  through  the  acoustic  axis.  Since  the  highest  frequency  of  interest 
was  500  hertz,  It  was  necessary  to  hold  a  positional  stability  to  within  approxi¬ 
mately  two  Inches  In  order  to  maintain  a  phase  stability  of  five  degrees.  Although  the 
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stiffness  of  the  boom  was  not  sufficient  to  guaranty  this  degree  of  stability  in  an 
environment  of  fluctuating  currents,  steady-state  measurements  revealed  no 
measurable  phase  fluctuations  which  could  be  attributed  to  hydrophone  motion. 
However,  attempts  to  repeat  measurements  after  the  boom  had  been  moved  resulted 
in  phase  differences  as  great  as  12  degrees,  part  of  which  might  be  attributable  to 
deviation  in  hydrophone  position,  with  the  remainder  probably  due  to  temperature 
variations.  The  resolution  of  the  phr«e-measuring  instrumentation  was  one  degree. 
The  range  from  the  center  of  the  array  to  the  hydrophones  varied  from  189  to  200 
feet,  depending  on  the  boom  angle.  These  ranges  provide  a  good  approximation 
to  far-field  conditions  for  the  magnitude  of  the  acoustic  field  in  the  major  lobe. 
Figure  10  illustrates  the  computed  400-hertz  directivity  patterns  in  a  vertical  plane 
perpendicular  to  the  plane  of  the  array  for  a  uniformly  vibrating  plane  rectangular 
piston  30  feet  wide  by  50  feet  high,  corresponding  to  the  dimensions  of  the 
ARTEMIS  source.  The  solid  lines  represent  the  computed  far-field  amplitude  and 
phase  patterns,  whereat  the  dashed  llnescfcplct  the  computed  near-field  patterns. 

The  far-field  values  were  computed  for  a  constant  range  and  referenced  to  zero 
decibels  in  magnitude  and  zero  degrees  In  phase  on  the  acoustic  axis.  The  near- 
fteld  data  were  computed  at  ranges  from  the  center  of  the  array  corresponding  to 
the  actual  hydrophone  ranges  and  then  referenced  to  the  same  range  at  the  far- 
field  patterns  In  order  to  provide  a  direct  comparison  between  the  values  which 
would  be  obtained  from  an  ideal  transducer  at  the  ranges  employed  and  in  the  far 
field.  It  can  be  seen  that  the  hydrophone  ranges  are  sufficient  to  approximate  the 
far-field  Intensity  to  within  0,5  decibel  for  angles  lest  than  ten  degrees  with 
respect  to  the  acoustic  axis.  However,  the  limited  range  does  distort  the  nulls 
considerably.  The  computed  near-field  phase  pattern  differs  from  the  far-field 
pattern  by  approximately  20  degrees  in  that  portion  of  the  major  lobe  located 
between  the  half  power  points.  This  phase  lag  results  from  the  Increased  range 
between  element  and  hydrophone  with  Increased  distance  of  the  element  from  the 
center  of  the  array  •  As  the  hydrophone  Is  moved  away  from  the  acoustic  axis, 
those  elements  having  decreasing  phase  lag  due  to  decreasing  range  contribute  an 
Increasing  share  of  the  intensity  and,  thus,  the  near-fleld  phase  exhibits  a  decreas¬ 
ing  lag  relative  to  the  far-field  phase.  In  Fig.  11,  the  difference  between  the 
computed  near-fleld  phase,  at  hydrophone  distance  on  the  acoustic  axis,  and  the 
equivalent  far-field  phase  Is  plotted  as  a  function  of  frequency  for  on  axis  and 
six  degrees  below  axis,  as  shown  by  the  solid  lines.  The  dotted  lines,  which 
represent  linear  phase  shifts,  are  plotted  for  comparative  purposes.  It  can  be  seen 
that  the  phase  error  caused  by  the  finite  hydrophone  range  is  very  nearly  a  linear 
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function  of  frequency,  representing  a  time  or  range  shift.  nnd  therefore  do**  r.st 
distort  the  form  of  the  transfer  function  In  the  vicinity  of  the  acoustic  axis.  It  can 
be  concluded  then  that  the  hydrophone  range  employed  was  sufficient  for  determin¬ 
ing  the  far-fleld  transfer  function  on  the  major  lobe.  It  is  also  apparent  that  precise 
measurements  for  angles  further  removed  from  the  acoustic  axis  require  a  greater 
hydrophone  range.  In  general,  the  monitor  hydrophone  should  be  placed  at  the 
greatest  practical  range  during  calibration  procedures.  The  widely  used  rule-of- 
thumb  range  criteria  of  ■fr/X,  where  t  is  the  maximum  dimension  of  the  transducer, 
is  not  sufficient  for  obtaining  precise  directivity  patterns.  It  is  evident  that,  in 
cases  where  the  depth  of  nulls  Is  of  Interest,  as  In  receiving  arrays  In  which  a  null 
is  steered  to  a  noise  spoke,  careful  attention  must  be  given  to  the  ranges  at  which 
measurements  are  made.  This  also  points  out  a  problem  in  null  steering  when  a 
noise  source  is  at  close  range. 

Having  satisfied  the  necessary  conditions  of  adequate  range  and  positional 
stability  for  the  monitor  hydrophone,  the  fidelity  of  an  acoustic  source  can  be 
determined  by  measurement  of  the  transfer  function.  It  Is  a  way  of  completely 
specifying  the  input-output  relations  for  a  system.  The  practical  Importance  of  the 
transfer  function  Is  that  It  can  be  determined  for  a  linear,  time-invariant  system  of 
any  degree  of  complexity  by  a  set  of  steady-state  measurements  using  sine  waves. 

The  techniques  and  equipment  needed  for  such  measurements  are  readily  available. 

There  Is  some  uncertainty  to  what  extent  high-power  acoustic  sources  can  be 
considered  linear  and  time-invariant  systems.  It  Is  known  that  such  systems  are  not 
strictly  linear  over  a  large  dynamic  range.  There  Is  reason  to  believe  they  are  not 
time-invariant.  If  a  transfer  function  for  such  a  source  Is  measured  using  sine  waves, 
this  transfer  function  may  not  determine  the  input-output  relations  for  all  types  of 
inputs.  For  this  reason,  It  was  decided  to  measure  the  transfer  function  of  the 
ARTEMIS  source  using  two  different  types  of  Inputs.  First,  the  transfer  function  was 
measured  using  sine-wave  inputs,  Second,  the  transfer  function  was  measured  by 
spectral  analysis  techniques  using  a  pseudorandom  signal  as  the  input.  A  comparison 
of  these  two  types  of  transfer-function  measurements  answers  In  part  the  question  of 
the  adequacy  of  sine-wave  measurements  for  the  ARTEMIS  source. 

Sine-Wave  Measurements  of  the  Transfer  Function  fer  the  ARTEMIS  Source 

The  ARTEMIS  source  Is  represented  by  the  block  diagram  (Fig.  12).  The  Input 
to  the  system  Is  the  input  voltage  to  the  console  of  the  power  amplifier.  The  output 
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of  the  source  Is  the  acoustic  pressure  measured  by  a  hydrophone  held  fixed  about 
190  feet  from  the  center  of  the  face  of  the  transducer  array.  The  transfer  fuiiviiuri 
H(uu)  1$  the  ratio  as  a  function  of  frequency  of  the  complex  phasors  P(uu)  and  V(u)  of 
the  acoustic  pressure  and  input  voltage,  respectively. 


In  the  transfer  function  H(uJ,  there  Is  a  factor  representing  a  phase  log  which  Is 
proportional  to  frequency.  This  corresponds  to  distortionless  transmission  due  to  a 
time  delay.  Part  of  this  phase  lag  Is  due  to  the  time  delay,  Tp  (time  required  for 
the  propagation  of  a  plane  wave  from  the  center  of  the  array  to  the  hydrophone). 
This  represents  a  factor  in  the  transfer  function  of 

“|U)Tp 

The  transfer  function  can  now  be  written  as 

m  «  h0(u)."|u,t°  ,  (2) 

where  Hq(u)  Is  interpreted  as  the  transfer  function  to  zero  range. 

The  method  of  measurement  of  H0(u)  is  Illustrated  In  Fig.  13.  An  oscillator 
produces  a  sine-wave  voltage  waveform  of  frequency  f,  which  Is  the  Input  to  the 
source.  A  zero-crossing  detector  detects  a  zero  crossing  of  this  voltage  which 
produces  a  pulse  which  Is  delayed  for  a  time,  tq,  and  then  starts  a  counter.  A 
zero  crossing  of  the  acoustic  pressure  Is  detected,  producing  a  pulse  which  stops 
the  counter.  The  counter  has  been  counting  a  frequency  360  times  f,  to  the  count 
It  the  phase  of  H0(u>)  with  a  resolution  of  ±1  degree.  The  amplitude  of  the  Input 
voltage  and  acoustic  pressure  was  measured  using  an  average  reading  vacuum-tube 
voltmeter.  Measurements  were  completed  at  five-hertz  Increments  of  frequency 
from  350  to  450  hertz . 

Figure  14  Illustrates  the  resulting  transfer-function  amplitude  and  phase  for 
the  ARTEMIS  source.  The  data  are  Illustrated  for  three  values  of  current:  8,  20, 
and  40  amperes.  For  any  one  current  value,  the  amplitude  displays  approximately 
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a  two-declbel  change  over  the  frequency  band,  whereas  the  phase  is  very  nearly 
linear  and  Is  very  nearly  free  of  distortion.  The  linear  phase  results  from  time 
shifts  In  the  amplifier  and  transducer  plus  any  propagation  delay  which  has  not 
been  completely  subtracted  out.  A  linear  phase  shift,  however,  does  not 
represent  distortion.  Considerable  nonlinearity  Is  evidenced  by  the  divergence 
of  the  data  between  8  and  80  amperes.  The  nonlinearity,  however,  was  not 
evidenced  In  spectral  measurements  In  which  the  harmonic  content  of  the  acoustic 
signal  was  found  to  be  more  than  40  decibels  below  the  fundamental  when  the 
array  was  driven  at  420  hertz.  Also,  It  will  be  shown  that,  when  bioad-band 
signals  were  used,  the  measured  transfer-function  phase  was  Identical  with  the 
phase  characteristics  measured  with  a  continuous-wave  Input  having  an  input 
current  equal  to  the  average  current  for  the  broadband  signal.  The  nonlinearity, 
therefore,  appears  to  be  a  function  of  average  rather  than  Instantaneous  power 
and  does  not  result  In  signal  distortion. 

The  amplitude  and  phase  characteristics  of  the  source  are  plotted  In  Fig. 

15  as  a  function  of  hydrophone  angle  relative  to  the  acoustic  axis.  For  comparison, 
the  computed  near-fleld  characteristics  are  shown  as  dotted  curves.  In  these 
plots,  the  amplitudes  and  the  computed  phase  have  been  referenced  to  zero  on 
the  acoustic  axis.  The  close  agreement  between  computed  and  measured  values 
Indicates  that  the  real  array  has  a  response  which  closely  resembles  an  Ideal# 
uniform-velocity  piston. 

Measurement  of  the  Transfer  Function  Using  Pseudorandom-Sequence  Signals 

Signal  Characteristics  -  A  maximal  length  linear  shift  register  sequence 
of  ones  anef zeros  was  used  to  modulate  a  400  hertz  square  wave.  An  output 
of  one  from  the  shift  register  does  not  change  the  carrier.  An  output  of  zero 
Is  equivalent  to  multiplying  the  carrier  by  -I  or  bhanglng  Its  phase  180  degrees. 

This  binary  waveform  was  then  passed  through  a  bandpass  filter,  United 
Transformer  Corp.  type  CF-765A,  having  a  center  frequency  of  400  hertz  and 
pats  band  from  345-455  hertz.  The  output  of  the  filter  was  the  signal  used  as 
the  Input  of  the  source.  This  signal  It  referred  to  at  phase  modulation,  although 
It  might  equally  well  be  called  suppressed-carrter  amplitude  modulation. 

The  signal  has  a  low  peak-to-rms  ratio. 

Two  pseudorandom-sequence  signals  were  used  In  the  transfer-function 
measurements.  Both  had  the  same  duration.  The  difference  was  the  shift-rate 
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and  feedback  connections  used  In  the  shift-register  generator.  Different  shift 
rates  were  used  to  change  the  bandwidth  of  the  signal .  The  power  spectrum 
of  the  signal  Is  roughly  a 


shape  centered  about  the  carrier  frequency .  The  first  nulls  In  the  spectrum  occur 
at  frequencies  of  the  carrier  plus  or  minus  the  shift  rate.  Table  3  gives  a 
description  of  the  two  signals. 


TABLE  3 


Carrier  frequency 

400  Hz 

400  Hz 

Duration 

41  sec 

41  sec 

Shift  rate 

100  Hz 

50  Hz 

Shift-register  length 

12  bits 

II  bits 

Frequency  of  first  nulls 

300,  500  Hz 

350,  450  Hr,  .  , 

Shift-register  sequence^ 

Degree  12 

Degree  II 

Code  M0I23F 

Code  147:52056 

6  Error  Correcting  Codes,  (W.  W,  Peterson) ^  the  Mil  Press,  IM 

Method  of  Measurement  -  The  cross-power  spectrum  between  two  signals 
having  Fourier  spectrums  X(uu )  and  Y (uj )  Is  defined  as 

Sx>)-X(u>)Y(u>)*,  (3) 

where  *  represents  complex  conjugation.  The  power  spectrum  of  X(ui)  Is 
defined  as 

Skx(ui)-X(u>)X(uu)*  «|X^)|2.  (4) 

The  spectrum  X(u))  Is  to  be  considered  as  the  input  to  a  linear  time-invariant 
system  and  Y(uu)  the  output  of  the  system.  The  system  has  a  transfer  function 
H(io),  so 

Y(ai)  ■  H(uj)X(uj),  (5) 
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and  substituting  this  Into  definition  of  cross-power  spectrum  gives 

S  to)  =X(u>)H(u.)*Xto)*  =  H(uj)*S  (tu).  (6) 

xy  xx 

From  Eq.  (6),  the  real  and  Imaginary  components  of  the  transfer  function  are 

(7) 


(8) 


To  determine  the  transfer  function  H{u>)  of  the  system,  It  is  necessary  to  measure 
the  real  and  Imaginary  components  of  the  cross-power  spectrum  between  input 
and  output  S^to)  and  the  power  spectrum  of  the  Input  S^to)* 

The  cross-power  and  power  spectrum  were  measured  by  the  analyzer  shown 
as  a  block  diagram  In  Fig .  16.  Each  block  represents  electronic  circuitry 
which  does  the  analog  of  such  operations  as  amplification,  multiplication, 
filtering,  adding,  and  Integrating .  A  measurement  of  cross-power  or  power : 
spectrum  using  this  analyzer  refers  to  an  average  over  a  narrow  band  of 
frequencies.  The  measured  cross-power  spectrum  about  frequency  u> „ 

(angular  frequency)  Is  denoted  by  S'  (ui  0)  and  is  deflnod  as 

*y 

%  +1  §  “b  +  I 

s'  to°)-|  f  B  s  toJcto-^J*  u  X(ui)Y  to)*cto  (9) 

where  B  is  the  analyzing  bandwidth,  vdth  the  prtme^indlcatlng  an  average  over 
the  band  of  frequencies  between  tu0  -  and  uo0  +  Similarly,  a  measurement 
of  power  spectrum  S'  (u^)  Is 

XX  Q 

+  § 

=  B  (,0) 

UJ0  -  2 


RetS  to)l 
Ro  [H  (uj  )]  =  m  j . 


-lm[S  (id)] 

ImtHto)]  =  T$n~ 
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Substituting  Into  Eq.  (9)  for  S 


xy 


(id)  from  Eq.  (6)  gives 


S'xyW  “  f  J*  H  (uj  )*Sxx(tu  )du> . 
«W-| 


(II) 


In  order  that  the  transfer-function  measurements  have  a  definite  mecnlng,  It 
must  be  assumed  H(ou)*  Is  constant  for  the  im  Interval  of  the  Integration  so  that 
we  can  write 


.  B 

U)0  +  2 

S'  s  to)*,. 


from  which 


S'x>ito.)  -  HtoJ*S 

The  transfer  function  Is  then  given  by 


Re[H  (uu  0)] 


(12) 


(13) 


(14) 


ImIHto ,)]  -  ~lmlS'xy <“'«>)  ,  (15) 

where  Re[S'  (a)  0)J#  lm[S*  (ai  0)] ,  and  S'  (u>  0)  are  the  quantities  to  be 

r  j.  *y  *y  * 

measured  by  the  analyzer. 

Equations  (9)  and  (10)  Indicate  spectrum  measurements  are  made  using 
bandpass  filters  of  bandwidth  B.  The  analyzer  shown  In  Fig.  16  actually 
uses  low-pass  filters  which  pass  frequency  components  between  -B/2  and 
+B/2.  It  will  now  be  shown  that  the  analyzer  using  low-pass  filters  measures 
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the  quantities  S'^(w  0)  and  S'^(uj  0).  The  gains  associated  with  the  various 

blocks  in  the  analyzer  can  be  ignored.  Since  the  transfer  function  is  a  ratio  of 
two  measurements  using  the  same  equipment,  the  gains  are  eliminated  from  the 
final  result.  It  Is  necessary  to  have  equal  gains  in  various  parallel  channels 
and  this  is  easily  obtained. 

The  Fourier  spectrum  of  the  Input  X(tu)  is  broken  up  into  a  positive  frequency 
part  X^(uu)  and  a  negative  frequency  part  X2(u>)* 

X(uu)  =  Xj(u>)  +X2(uu),  (16) 

where 

Xj(uj)  =  X(uu)  uj>-0 
X,(u>)  =  0  to  <0 

and 

X2(cu)  “0  uj  >0 
X2(uj)  3  X(uu)  ui<0. 


A  similar  definition  is  used  for  V(ui).  Referring  to  Fig.  16,  the  input  time 
function  x(t),  corresponding  to  X(w),  is  multiplied  by 

-K* 

COSU)  0t  *»  - J - 

In  the  block  labeled  modulator.  The  resultant  time  function  Is  passed  through 
a  low-pass  filter  which  passes  angular  f  •  /ency  component  within  ±13/2  of 
zero  frequency.  Using  the  relation  for  t  'plex  modulation  of  a  time  function, 
the  spectrum  of  the  output  of  the  low-pass  filter  Xn(uu)  Is 

X,(<D-U)0)  +  X. (uu+uu  0) 

Xp(uu)=-i - j - ! -  for 

X  (ou)  =  0  elsewhere. 

P 


<U) 


<1 


and 


07) 
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The  input  x(t)  is  also  multiplied  by 
ei“>  oLb’1»  o* 


sinuu  J  » ■ 


TT 


and  low-pass  filtered,  resulting  In  a  waveform  with  spectrum  X  (uu). 

q 


X  (iu'-uu0)-X,(uj  +  uj0) 

*,«•>  .  -i - ^ - 

X  (uu )  =  0  elsewhere. 

q 


B  B 

for  -  -j<u»  and 


08) 


Treating  the  spectrum  Y(uu)  the  same  way  gives 

Y-(u)  -  u>  0)  +  Y ,  (u)  +  («  0)  q  » 

Y_(uu )  -  —  . . .  for  -  j-  <uj  <  j  ,  and 


P-  •  r 

Y2(u>  -  ud0)  -Y^id  +tt)0) 


V“> 


for  -  <jjj  . 


(19) 

(20) 


Solving  Eqs.  (17)and  (18)  forX^(u)  +uj0)  and  (19)  and  (20)  for  Yjtw  +  u>0) 
gives 

X | (ui  +  0)  "  Xp(ou)  -  |X^(w)  (21) 

Y^+ui,)  -  Yp(w)  -  IYq(uu)  (22) 

B  B 
for  -  7j|<UU  <TJf  » 

Returning  to  the  definition  of  measured  cross-power  spectrum  (Eq.  9). 

J  B 

.  +  5 

S'x  (uip)  *  g-  J*  g  X (uu )Y (uu )*duu ,  and,  since  the  Integration  Is  over 

X/  8  1 
positive  frequencies, 
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,  B 

<".+■5 

S'  «»J  -i  J  B  X^Tf^rdm. 

y  ».-* 


/*o\ 


Making  change  of  variable*  give* 


X^U)  +  u)  0)Y ,  (id  +  tD0)*dtD. 


(24) 


Substituting  Into  Eq.  (24),  X^  +  iu#)  and  Y,  (u>  +  id,)*  from  Eq».  (21)  and 
(22)  gives 


frO  "  J* 


B/2 


B/2 


[X  (id)  -  |X  (u>)!  tY >)  -  (Y  (<D)]*d<D . 
P  *1  P  “ 


(25) 


Expanding  and  rearranging  term*,  thli  become* 


Bj/2 

+1  f  [X  (id)Y  (id)*  -X  (id)Y  >)*]diD.  (26) 

-  B/2  p  q  q  P 

The  value*  of  both  definite  Integral*  In  Eq.  (26)  are  real.  Thl*  follow*  from 
the  fact 

X  (id)  -  X >  u>)*,  (27) 

P  P 

X  (id)  -  X  (-id)*, 

9  ^ 

X  (id)  -  Y  (-  id)*,  and 

P  P 

Y>)  -  Y  (-»)*, 

*1  " 
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which  can  be  shown  by  manipulation  or  noting  that  they  represent  Fourier 

spectrums  of  real  time  functions.  Therefore,  the  first  integral  in  Eq .  (26) 

represents  Re[S'  (uu0)l  and  the  second  Integral  ImlS1  (a>  0)1 .  The  analyzer 
xy  xy 

performs  the  Integration  over  frequencies  Indicated  by  Eq.  (26)  by  multiplying 

the  corresponding  time  functions  and  Integrating.  If  fj(t)  and  fj(t)  are  real 

time  functions  and  F ^  (uu )  and  Fj^)  ore  their  Fourier  transforms, 


00  00 

J*  U0Ut)dt  =  J  F1(ut))F2(uu)*duu. 
■■00  —  00 


(28) 


The  final  value  of  the  output  of  the  Integrators  Is,  respectively,  Re[S'  (uu  0)] 

xy 

and  lm[S*  (uu  J] .  To  measure  the  power  spectrum  of  X(tu)  using  the  same 
xy 

analyzer,  putY(u>)  «=X(tu). 


Measurement  of  the  transfer  function  of  the  ARTEMIS  array  using  the  cross¬ 
power  spectrum  analyzer  is  shown  in  Fig.  17.  The  generator  output  Is  bandpass 
filtered  to  345-455  hertz,  producing  a  signal  with  spectrum  V(tu)  used  as  the 
Input  to  the  source.  The  generator  output  Is  precisely  delayed  by  a  time  t 

— |(U  T  ® 

and  bandpass  filtered,  producing  a  signal  with  spectrum  V(w)e  1  o.  This 
signal  is  used  as  the  reference  input  to  the  analyzer  or  what  was  called  X(u>) 

In  the  previous  discussion.  The  acoustic  pressure  at  the  hydrophone  P(uu )  Is 
the  Y  (ai)  Input  to  the  analyzer.  The  analyzer  measures  the  transfer  function 
of  the  source  H 0 (a)  0)  as  defined  In  Eq.  (4).  H,(u)0)  is  the  transfer  function  of 

the  source  at  angular  frequency  U>  0  after  the  factor  e"^  °T  due  to  pure  time 

delay  has  been  eliminated.  The  elimination  of  e"^°T  factor  by  Introducing 
a  delay  reference  Into  the  analyzer  ts  an  essential  part  of  the  measurement. 

The  delay  causes  a  phase  shift  of  about  14  degrees  per  hertz  frequency 

change.  If  this  phase  shift  Is  not  eliminated,  It  Is  not  possible  to  assume  the 
measured  transfer  function  Is  constant  over  the  analyzing  bandwidth.  The 
analyzing  bandwidth  used  In  the  measurements  was  four  hertz.  Measurements 
were  made  at  frequencies  tour  hertz  apart  from  350  to  450  hertz. 
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Another  type  of  measurement  was  made  to  cheek  the  coherence  of  the 
ARTEMIS  source.  Because  of  the  large  time-bandwidth  product  of  the  i., 

—  *  ^ _ _ _  i. .  / .  I.  _  .  .s.  AAA\  il.  ......  •  .  • 

ni«uaui«fiimiii»  \guuMi  auv/^  iiio  ctu»»~puwai  epocuuin  m«a»uromonr  ignore* 

Incoherent  energy.  By  measuring  both  the  Input  and  output  power  spectrum 

S1  (te)  and  S'  (u>),  It  Is  possible  to  define  a  coherence  ratio  for  the  source, 
xx  yy 


Y(w  c) 


(29) 


sVy<“'.) 


This  ratio  Is  equivalent  to  the  normalized  cross-carreltitlan'fbr.ithe  band  of 
frequencies 


»e 


B 

“2 


<  tl)  <U) 


e 


+ 


B 

T 


If  the  system  Is  fully  coherent,  Y(u>«)  should  equal  unity  at  all  frequencies. 
Values  less  than  unity  are  Interpreted  In  terms  of  loss  of  coherence. 


Results  -  The  transfer-function  phase  data  for  B  and  60  amperes  sine  wave 
shown  in  Pig.  14  Is  reproduced  to  an  expanded  scale  In  Fig.  18,  along  with 
the  phase  data  measured  with  a  pseudorandom  sequence.  The  average  value 
of  the  rms  current  for  the  pseudorandom  sequence  signal  was  80  amperes. 

The  signal  had  a  shift  rate  of  100  hertz.  The  excellent  agreement  between 
the  broadband  and  continuous-wave  data  at  80  amperes  Indicates  that  the 
phase  nonlinearity  of  the  source  Is  a  function  of  average  rather  than  Instantaneous 
current  since  the  Individual  frequency  components  of  the  pseudorandom 
sequence  signal  have  amplitudes  whteh  are  very  small  compared  b  80  amperes. 

It  would  not  be  expected  then  that  the  nonlinearity  of  the  source  would  degrade 
Its  performance  with  broadband  signals. 

The  magnitude  and  phase  characteristics  of  the  transfer  function  at  angles 
of  0s,  7.5°,  and  20°  relative  b  the  acoustic  axiri  are  plotted  In  Figs.  19 
and  20.  The  7.5*  position  Is  below  the  acoustic  axis  on  one  side  of  the  major 
lobe  and  the  20*  position  Is  below  the  acoustic  axis  on  the  first  minor  lobe. 

The  solid  lines  were  plotted  from  data  of  the  100  hertz  shift-rate  sequence 
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and  the  dotted  curves  with  circled  points  from  the  50  hertz  shift-rate 
sequence.  The  phase  data  has  been  referenced  to  zero  degrees  at  400  hertz 
by  subtracting  a  linear  phase  shift.  The  phase  and  amplitude  distortion  Is 
small  at  all  artgies  investigated,  and  no  significant  difference  !» observed 
between  the  two  types  of  modulated  waveforms. 

The  coherence  ratios  are  not  Illustrated  since  they  have  no  significant 
deviation  from  unity.  For  all  frequencies  within  the  operating  band  and  for 
all  angles  with  respect  to  the  acoustic  axis  from  0  to  -20°,  the  coherence 
ratios  computed  from  the  measured  data  fell  within  the  range  of  0.99  to 
1*01 .  Since  the  accuracy  of  the  spectral  analysis  equipment  Is  considered 
to  be  ±2  percent,  the  coherence  ratios  are  unity  within  the  accuracy  of 
measurement. 


Cross-Correlation  Function  Measurements 

Figure  21  shows  the  method  of  measurement  of  the  cross-correlation 
function  between  the  Input  and  output  of  the  source.  The  Input  was  a 
pseudorandom-sequence  signal  generated  by  the  digital  generator  previously 
described.  The  output  of  the  digital  generator  was  sampled  at  a  rate  of 
800  times  a  second  synchronous  with  tts  shift  pulses.  These  samples  were 
entered  Into  one  end  of  a  64-bit  shift  register.  By  picking  off  the  signal 
along  the  stages  of  the  shift  register,  a  time  delay  t,  In  Increments  of 
1.25  milliseconds,  from  0  to  80  milliseconds  was  obtained.  By  varying  the 
sampling  of  the  output  of  the  generator  relative  to  Its  clock  from  0  to  1 .25 
milliseconds,  any  delay  from  0  to  80  milliseconds  can  be  obtained.  The 
output  of  the  generator  and  shift  register  were  filtered  by  a  matched  pair  of 
bandpass  filters  (United  Transformer  Corp.  type  CF-765A)  to  obtain  the 
signal  v(t)  Into  the  source  and  a  reference  signal  v(M)  used  In  the  cross- 
correlation. 

As  shown  In  Fig.  21,  the  acoustic  pressure  p(t)  was  multiplied  by  v(t-i) 
and  Integrated  for  the  41-second  duration  of  the  signal  •  The  final  value  of 
the  Integrator  was  held  and  read  out  using  a  digital  voltmeter. 

Aside  from  a  constant  factor  which  represents  the  various  gains  In  the 
measurement  system,  the  final  value  of  the  integrator  output  Is  the  cross¬ 
correlation  R(t)  between  the  input  and  output  of  the  source. 

T 

R(t)  -  J  p(t)v(t-i)dt.  (30) 
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To  reduce  the  time  required  to  complete  the  measurements.  It  was  decided 
to  use  only  values  of  t  which  correspond  to  the  peaks  (positive  ar  negative) 
of  the  fine  structure  of  the  cross-correlation  function.  Peaks  occur  for  values 
of  t  1 ,25  milliseconds  apart;  that  is,  half  the  period  of  the  400  hertz  carrier. 
By  varying  the  sampling  of  the  generator  output  over  1 .25  milliseconds 
relative  to  the  carrier,  a  few  measurements  would  locate  a  peak  In  the  fine 
structure.  After  that  was  found,  other  peaks  could  be  obtained  by  tapping  off 
the  signal  along  the  shift  register. 

Since  there  was  a  pure  time  delay  of  about  40  milliseconds  In  the  source 
and  a  delay  of  up  to  80  milliseconds  was  available,  It  was  possible  to  obtain 
values  of  the  cross-correlation  function  40  milliseconds  on  either  side  of  the 
peak  of  Its  envelope.  This  range  of  t  determines  the  cross-correlation  through 
Its  main  lobe  out  past  the  first  side  lobes. 

The  normalized  cross-correlation  function  Is  defined  as 


o(t) 


T 


,r  p(t)v(t-T)dt , 

o 


(31) 


2  2 

The  Integrals  of  p(t)  and  v(t)  were  obtained  by  the  same  multiplier-integrator 
combination  used  In  the  cross-correlation  function  measurements.  It  Is  then 
possible  to  compute  o(t). 


The  correlograms  obtained  between  the  Input  signal  and  the  hydrophone 
outputs  confirm  that  negligible  decorrelation  Is  produced  by  the  ARTEMIS 
source.  Figure  22  illustrates  the  normalized  cross-correlations  obtained 
with  50-  and  100-hertz  shift-rate  sequences.  There  ere  no  measurable 
differences  between  these  curves  and  plots  of  the  normalized  autocorrelation 
functions  of  the  corresponding  Input  signals.  Although  these  plots  were  made 
from  data  obtained  with  a  hydrophone  on  the  acoustic  axis,  they  are  identical 
with  correlations  obtained  at  angles  up  to  22  degrees  below  the  axis.  It  can 
be  concluded  that  this  source  generates  an  acoustic  f| «ld  which  Is  a  replica 
of  the  Input  time  function  with  negligible  distortion. 
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SUMMARY 

The  potential  gain  In  sonar  performance  to  be  obtained  by  maximizing 
signal  power  justifies  Increased  effort  addressed  to  this  problem.  Adequate 
techniques  are  needed  to 

1 .  Determine  the  factors  which  limit  the  allowable  Input  power  to 
acoustic  sources, 

2.  Test  nondestructive^  the  power  limitation  In  existing  transducers 
when  driven  with  useful  types  of  modulated  signals, 

3.  Establish  design  criteria  for  optimizing  power  capability  of  sources 
to  be  used  with  specific  types  of  signals,  and 

4.  Test  transducers  for  phase  and  amplitude  distortion. 

Acoustic  sources  having  dimensions  which  are  large  compared  to  the  sound 
wavelength  In  water  are  Inherently  subject  to  large  variations  In  acoustic 
loading  over  the  radiating  face.  Variations  In  loading  resulting  from  acoustic 
interactions  are  both  spatial  and  temporal,  being  related  to  the  physical 
geometry  of  the  source  and  the  time  function  of  the  signal  modulation.  If 
the  range  of  values  over  which  the  real  and  reactive  parts  of  the  acoustic 
load  can  vary  is  determined  fora  particular  array,  and,  (fall  elements  of 
the  array  can  be  assumed  to  be  electrically  and  mechanically  Identical,  then 
the  power  capability  of  the  array  can  be  predicted  from  an  analysis  of  a 
single  element.  The  transient  response  of  the  power-1  Imltlng  factors  within 
a  transducer  element  are  difficult  to  predict  analytically  or  measure 
experimentally.  These  factors,  which  might  be  In  the  form  of  mechanical 
or  electrical  stress,  temperature,  or  dimensional  tolerances,  are  usually 
peak  limited.  Each  factor  has  an  upper  limit  which,  If  exceeded,  results 
In  nonlinear  operation  or  failure  of  the  element.  Which  factor,  however, 
actually  limits  the  power  input  to  the  element  depends  upon  external  conditions 
such  as  acoustic  load,  excitation  frequency,  or  driver  output  Impedance. 

The  limiting  factor  under  one  set  of  conditions  may  be  replaced  by  another 
factor  under  a  different  set  of  conditions. 

Analog  modeling  Is  an  attractive  technique  for  transient  analysis  of  the 
power-limiting  factors  In  that  all  dynamic  quantities  ere  represented  as 
voltages  or  currents  and  thus  are  readily  accessible  for  measurement  and 
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external  variables  such  as  acoustic  load  are  easily  controllable.  This 
method  Is  limited,  however,  by  the  ability  to  faithfully  represent  a 
mechanical  system  with  all  Its  decrees  of  freedom  by  an  electrical 
circuit.  In  this  method,  as  In  any  other.  It  is  necessary  to  be  able  to 
predict  the  range  over  which  the  components  of  the  acoustic  load  can  vary 
in  operation.  No  adequate  method  Is  presently  available  to  do  this. 

The  Increased  emphasis  on  coherent  signal -processing  techniques  has 
brought  about  a  need  for  more  sophisticated  transducer  test  procedures  than 
have  been  conventionally  required.  Techniques  used  In  obtalnii\|  the  transfer 
function,  coherence  ratios,  and  cross-correlation  between  the  electrical 
Input  and  acoustic  output  of  the  ARTEMIS  acoustic  source  serve  as  an 
example  of  the  types  of  Instrumentation  required. 
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Fig.  2  •  Basic  configuration  of  an 
illustrative  ceramic  element 


L0  ■  ACTUAL  SHUNT  TUNING  INDUCTANCE  IN  HENRIES 

Co  ■  ACTUAL  SHUNT  CAPACITANCE  OF  CERAMIC  STACK  IN  FARADS 

Ci  «  EQUIVALENT  CAPACITANCE  DUE  TO  CERAMIC  STIFFNESS 

C8  •  ••  "  «  •'  PRESSURE  RELEASE  MAT'L 

L,  ■  EQUIVALENT  INDUCTANCE  DUE  TO  MASS  mi 

L|  ■  "  •'  ••  »  •>  mi 
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Ci  «  *  ■  ELECTROMECHANICAL  TURNS  RATIO 

Ca  ■  flVia  .  , 

ii  S  Mz  ARE  STIFFNESS  COEFFICIENTS  FOR 
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Fig.  3  -  Equivalent  circuit  of  the  element 
shown  in  Figure  2 
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Z.  ■  EQUIVALENT  IMPEDANCE  DUE  TO  THE  ACOUSTIC 
RADIATION  IMPEDANCE 

Cmi  *  EQUIVALENT  CAPACITANCE  DUE  TO  OUTER  MASS  mi 
Cm2  *  EQUIVALENT  CAPACITANCE  DUE  TO  INNER  MASS  m2 
Rm  *  EQUIVALENT  DISSIPATION  ELEMENT 
Lm  *  EQUIVALENT  INDUCTANCE  DUE  TO  SPRING  STIFFNESS 
R,  ■  ACTUAL  WINDING  RESISTANCE 
L|  *  ACTUAL  WINDING  INDUCTANCE 
C,  ■  SIGNAL  COUPLING  CAPACITOR 
R  ■  SOURCE  RESISTANCE 

Fig.  8  -  Equivalent  circuit  of  the  Astemia  transducer  element 
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imam  allowable  power,  input  to  the  Artemis  source  as  determined 
experimentally  and  by  analogue  circuit  analysis 


PHASE  -  DEGREES  INTENSITY 


PUTED  FAR  FIELD  PATTERN 
PUT  ED  PATTERN  AT  HYDROPHONE 


0  2  4  6  8  10  12  14  16  IB  20  22  24 

HYDROPHONE  DEPRESSION  ANGLE  WITH  RESPECT  TO  ACOUSTIC  AXIS- 

DEGREES 

Fig«  10  -  Computed  400  hertz  directivity  pattern  for  a  piston 
having  the  dimensions  of  the  Artemis  source 
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Fig.  12  -  Block  diagram  of  the  Artemis  source 


Fig,  14  -  Transfer  function  of  the  Artemie  source  for  three 
value*  of  aine  wave  current 
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Fig.  15  -  Acouatic  amplitude  and  phase  characteristics 
of  the  Artemis  source 
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FREQUENCY-Hz 

Fig.  18  -  Comparison  of  transfer  function  phase  of  the  Artemis 
source  as  measured  with  sinusoidal  and  pseudo-random  signals 
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Fig.  20  -  Phase  of  the  Artemis  source  transfer  function  for  three 
angles  relative  to  the  acoustic  axis 


Fig.  21  -  Instrumentation  used  for  measurement  a£  the  cross  correlation 
function  between  the  input  and  output  of  die  Artemis  source 
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The  generation  of  a  high- intensity  acoustic  field  which  is  a  faithful  replica,  or  at 
least  a  known  function,  of  the  input  signal  is  a  prime  requirement  in  modern  active 
sonar  systems.  Most  sonar  signals  in  use  are  modulated  sine  waves.  Therefore, 
transducers  must  be  capable  of  efficiently  transforming  this  type  of  signal.  It  is  evi¬ 
dent  that  design  criteria  and  teat  procedures  baaed  on  steady-state  sine-wave  signals 
are  not  adequate.  The  potential  gain  in  sonar  performance  to  be  realised  by  maximis¬ 
ing  signal  power  without  distortion  merits  increased  effort  addressed  to  this  problem. 

This  report  is  intended  to  point  out  some  of  the  problems  associated  with  the 
generation  of  maximum-energy  acoustic  fields  by  large  multielement  sources  and  to 
euggeat  techniques  of  analysis  and  testing.  The  use  of  electrical  analog  analysis  of 
power-limiting  factors  in  transducer  elements  appears  to  be  an  attractive  method  if 
a  means  can  be  developed  for  estimating  the  range  of  acoustic-loading  values  experi¬ 
enced  by  elements  in  an  active  array. 

A  description  of  calibration  procedures  used  and  results  obtained  in  tests  of  the 
ARTEMIS  acoustic  source  are  presented  to  serve  as  an  example  of  proposed  test 
methods , 
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